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04 ' The excited states of ^^Be, ^"^Be and ^^B have been studied with a method of antisym- 

}_( ' metrized molecular dynamics. In the predicted excited states we find novel molecule-like 

I structures with very exotic clusters such as ®He+*He in ^^Be and ^He+^Li in ^^B. The 

• origin of the ®He cluster development in the very neutron-rich nuclei is understood by the 

' new-type correlation among 4 neutrons and 2 protons. In this paper we also present our 

in I recent challenge to study sd-shell nuclei. Shape coexistence problems in '^'^Ar and '"^Ca are 
discussed. 

^ ' §1. Introduction 

o ■ 

I In the light nuclear region, a clustering aspect is one of the essential features in 

■T^lj- ' unstable nuclei as well as in stable nuclei. Owing to the progress of the experimental 

. techniques, information for the excited states of light unstable nuclei have been 

I increased rapidly. The exotic clustering of light unstable nuclei is one of the attractive 

subjects in the experimental and theoretical researchesL 

For example, theoretical studies on ^*'BeEP'l3)'3)'EP'El) suggested molecule-like 
structures in the excited bands K = 1^ and K = O2 ■ The other candidates for 
^ ' the molecule-like states are the excited states of ^^Be which have been recently dis- 

Ch ■ covered in the break-up reactions into ^He-|-^He and ^He-|-^HeQ). On the other hand, 

the structure of the ground state of ^^Be has been interesting because there are some 
experimental data which indicate the vanishing of magic number 8 in ^^BelZP'S*. How- 
ever ^^Be has not been studied enough by full microscopic calculations assuming no 
\ cluster core. It is also important to search for molecule-like states in other neutron- 

rich nuclei in order to understand universal features of clustering in unstable nuclei. 
Our first aim is to study the states covering the ground and the excited ones of light 
neutron-rich nuclei with antisymmetrized molecular dynamics (AMD). We study the 
structures of ^^Be and challenge to discover possible molecule-like states with exotic 
clusters in heavier neutron-rich nuclei, ^"^Be and ^^B. 

We apply the microscopic method of antisymmetrized molecular dynamics (AMD) 
which has already nroved to be a very useful approach for the structure study of 
general light nuclei lP In the study of ^■^cS', the author has proposed 

a new version of the method, variation after spin-parity projection (VAP) in the 
AMD framework which is very useful to investigate excited states with various kinds 
of structures such as spherical shell-model-like structures and clustering structures. 
The author and her collaborators have succeeded to describe the structures of excited 
states of ^"BeEl^ with the VAP calculations in the framework of AMD. AMD is the 



X 



typeset using pTpJ^sty <ver.0.8> 



2 



method that is very suitable to search for the possible exotic clusters in foreign nuclei 
because we do not need any model assumptions such as the existence of clusters. 

As mentioned above, clustering is an important feature in very light nuclei. 
However it is an open problem whether or not clustering effects are seen also in 
unstable sd-shell nuclei. I challenge to study structures of sd nuclei with AMD 
method. In this paper, the recent hot subjects of shape coexistence problems in 
^^Ar and ^'^Ca have been studied with AMD. Although they are stable ones, the 
studies of stable nuclei are very essential to investigate the unstable sd-shell nuclei. 

In this paper, the structures of the excited states of light neutron-rich nuclei, 
^^Be, ^^Be and ^^B are studied in Sec.^. We discuss the mechanism of the clustering 
development from the view point of correlation among nucleons in the single particle 
orbits. In the study of heavier nuclei with AMD, the shape coexistence problems in 
^^Ar and ^'^Ca are investigated in Sec.^ 

§2. Formulation 

An AMD wave function is a Slater determinant of Gaussian wave packets; 
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where the centers of Gaussians Zj's are complex variational parameters. Xi is an 
intrinsic spin function represented by which is also a variational parameter speci- 
fying the direction of the i-ih. intrinsic spin, is an isospin function which is fixed 
to be up(proton) or down(neutron) in the present calculations. 

In order to obtain wave function of an excited states of light nucleus, we vary 
the parameters Zj and = 1 ^ A) to minimize the energy expectation value for 
the parity and total angular momentum eigenstate (VAP calculations). 
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where the operator of total angular momentum projection P^jj^/ is / dODfjj^, {f2)R{f2). 
The integrationfor Euler angle i? is calculated numerically. We adopt the frictional 
cooling methodB^ to obtain the minimum energy states. Thus we can obtain the low- 
est state for a given spin parity with VAP calculations. For higher excited states 
we perform the energy variation for the orthogonal component to the lower states by 
superposing wave functions. More details of the AMD method for the excited states 
with the variation after spin-parity projection are described in Refsll3)'EP. By making 
VAP calculations for the lowest and the higher excited states with various sets of 
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total spin and parity {J^}, we obtain a lot of AMD wave functions {^^i, • • • ,^m}, 
which approximately describe the intrinsic states of the corresponding states. 
The number m is the number of the considered states. Final results are attained by 
diagonalizing a Hamiltonian matrix {P^^,<Pi\H\P^,^„<Pj) = 1 ~ m). 

For the structure study of heavier nuclei such as '^^Ar and ^'^Ca, we make vari- 
ation after parity projection but no spin projection (VBP) with a constraint AMD 
instead of full VAP calculations to save computational time. The adopted constraint 
in the present calculations is that the expectation values of total oscillator quanta 
must equal to a given number as (aa^) = W. After VBP calculations with a con- 
straint, we make spin projection to obtain a energy curve as a function of total 
oscillator quanta. The obtained states are superposed by diagonalizing Hamiltonian 
Matrix (P/^=^,^i|i/|P/^^„#,). 

§3. Results of i^Be, ^^Be and ^^B 

We apply the AMD method for the excited states of ^^Be, ^''Be and ^^B. The 
adopted interactions in this work are the central force of the modified Voikov No.l 



with case 3c3^ , the spin-orbit force of GSRSll^* and the Coulomb force. The Majorana 
parameter used here is m = 0.65, and the strength of GSRS force is chosen to be 
ui = —U2 = 3700 MeV. We choose a optimum width parameter i/ for the Gaussians 
of the single particle wave functions of each nucleus. With these parameters, the 
parity inversion of the ^^Be ground state can be reproduced. The binding energies 
of i^Be, i^Be and ^^B are 61.9 MeV, 59.7 MeV and 73.1 MeV, which are smaller 
than the experimental data 68.65 MeV, 69.77 MeV, 88.19 MeV, respectively. We 
have checked that the underestimation with the present interaction parameters can 
be improved easily by using smaller Majorana parameter as m = 0.61. With this 
parameter set the binding energy of ^^Be is 67.7 MeV which weh agrees to the 
experimental data. We also find that the excitation energy of 0^ in ^^Be calculated 
with m = 0, 61 is 4.2 MeV which is almost as same as the one with m = 0.65. 

In the calculated results of ^^Be with AMD, a lot of excited states appear in 
the low-energy region. The energy levels of ^^Be are presented in Fig. 1. The 
theoretical levels of 43", 6^ and 8^ states well correspond to the recently observed 
excited states H'. By analyzing the intrinsic AMD wave functions, we can classify the 
excited states into rotational bands such as K'" = 0^, 02^, O3", 1^". The interesting 



point is that the newly observed levelsEF at the energy region above 10 MeV have 
been found to belong to the third rotational band = 03". It is also interesting to 
see the ground = ^'^ band properties because the vanishing of the neutron magic 
number N = 8 occurs in ^'^Be. Even though the ^'^Be is a neutron magic nucleus 
with N = 8, the intrinsic state of the ground 0"^ state is not the ordinary state 
with the closed neutron p-shell, but a prolately deformed state with a developed 
clustering structure, which is dominated by 2p — 2h configurations in terms of single- 
particle orbits. As a result, the ground K'^ = band starts from the ground 0"*" 
state and reaches the band terminal at J'^ = 8"*" with the highest spin state in the 
2huj configurations. On the other hand, the neutron p-shell closed states with Ohuj 
configurations construct the second K'" = 0^ band which consists of the second 
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0^ and the second 2^ states. For an experimental evidence, the strength of the /? 
decay from ^^Be(O^) to ^^B(l^) is helpful to estimate the breaking of the neutron 
p-shellQ\ As is expected, the theoretical value of B(GT)=0.8 is enough small as the 
experimental data B(GT)=0.59 because the component of the p-shell broken state 
in the parent ^^Be(O^) makes the expectation value of Gamow- Teller operator to 
be small. It is consistent with the pioneer works EJ)' LP. It should be pointed out 
that this is the first theoretical work which can reproduce systematically the energy 
levels from the ground state to the highly excited states. We have predicted the new 
K'^ = OJj" rotational band from 0^, 2^, 62 ^^'^ ^1 states. The state has 
the mostly developed ^He+^He clustering. The developed clustering becomes weak 
gradually with the increase of the total-angular momenta J, and finally it changes to 
the same spin aligned state at = 8^ as the one in the K = Of band. The reason 
for the spin alignment at = 8^ is because the states in the = O'^ band are 
the other 2huj states which must have the same highest spin state as the one in the 
ground band at the band terminal 8^. Although the other evidence for thejpagic 
number vanishing is the low-lying 1~ state which has been measured recently eP, the 
excitation energy of 1~ state is overestimated in present calculation. The mechanism 
of energy gain of negative parity states has been still an open problem. 



> 20 



(U 

W 10 
c 

X 

W 



(8+) 


+—8+ 




(6+) 






(4+) 


3—2+ 


_& 

4' 




4+_ 


2+r=" 

3" 





2+' 

0+ 


0+ 


EXP 


AMD 



Fig. 1. Excitation energies of the levels in ^^Be. 
Theoretical results are calculated by VAP cal- 
culations based on AMD. Experimental data 
iken from the Table of Isotopes and Refs 
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Fig. 2. Excitation energies of the levels in ^*Be 
calculated by VAP calculations based on 
AMD. Excitation energies are plotted as a 
function of J( J + 1) where J is a total spin of 
the state. 



By analyzing the single particle wave functions in the intrinsic state, the mech- 
anism of the clustering development in the first 0"*" state can be understood by an 
idea of molecular orbits surrounding 2n! cores. Although the molecular orbits in ^^Be 
have been suggested theoretically li^'Q^, the important point of the present work is 
that the structure of 2a cores and surrounding neutrons in ^^Be is formed auto- 
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matically in the energy variation though the existence of any clusters nor molecular 
orbits is not assumed. On the other hand, it is natural to consider the third O"*" state 
as ^He+^He clustering structure instead of the a and surrounding neutrons because 
the distance between two clusters is too far to be described in terms of the molecular 
orbits. In this case, an a cluster goes outward far enough to form a ^He cluster in 
the correlation with 2 valence neutrons. 

Next we present the results of ^^Be. Figure2 shows the energy levels of the posi- 
tive parity states of ^^Be. While few excited states are known experimentally, many 
excited states are predicted in the theoretical results. There are very recent experi- 
mental datallZl of a few levels just above the threshold energy(9.1 MeV excitation) 
for the separation into ^He and ^He. 

The excited states in the region J < 6 are classified into 3 rotational bands, 
X = 0^, = 0^ and X = 2"*". In most of the states in ^^Be, 2a cores are formed in 
the calculated results. As a result, 0^ and states have prolate deformations with 
the deformation parameters /3 = 0.49 and /3 = 0.64 respectively. By analyzing the 
single particle wave functions of the intrinsic states, we found that all the states in the 
K = and K = 1'^ bands are dominated by the normal 0?ia; configurations. What 
is interesting is the exotic clustering structure of the excited states in the second 
K'^ = 0^ band which comes from 2fiu} configurations. In the present calculations, 
this band is predicted to start from the second = 02^ state at about 5 MeV and 
to reach the 8^ state. They have the well-developed ^He and ^He clusters which are 
very neutron-rich nuclei themselves. The surface cut for the matter density p > 0.16 
nucleons/fm^ of the intrinsic state of 02^ state is shown in Fig. 3(a). The ^IIe(^He) 
cluster is clearly seen in the right (left)-handed side in Fig. 3(a) where the surrounding 
low density region is omitted. As seen in the figure, the spatial clustering develops 
remarkably. It is interesting that we can find an origin of the clustering in the single 
particle wave functions. According to the analysis of the energies and behaviors 
of the single particle wave functions, the highest 4 neutron orbits correspond to 
sd-like orbits in the deformed system. Roughly speaking, a pair of spin up and 
down neutrons occupies a higher spatial sd-like orbit and the other pair of neutrons 
occupies the lower spatial orbit. Fig.3(b) and Fig.3(c) shows the two types of the 
spatial orbits for the 4 valence neutrons which contain more than 80 % components 
of positive parity states. It is found that the 4 neutrons occupy the sd-like orbits 
modified in the prolately deformed system. If we call the longitudinal direction of the 
prolate deformation as z-axis, the highest nucleon orbit seen in Fig. 3(b) is associated 
with the sd orbit with a form of yzexp[— z^r^] which has nodes for the rotation at 
the X-axis. On the other hands, the spatial orbit for the lower one is similar to the 
orbit of exp[— i/r^] which has nodes along (Fig.3(c)). These neutron orbits 

are stabilized due to the development clustering structures because they gain their 
kinetic energies of the nodes along the z-axis. From such a viewpoint of the single 
neutron orbits we can consider that four neutrons in sd-like orbits and 2 protons in 
p-like orbits forms a ^He cluster out of a '^He core. In another word, the extremely 
developed ^He cluster is originated from a new-type correlation among 4 neutrons 
in sd shell and 2 protons in p shell. 

In the theoretical results of ^^B, we find an exotic clustering ^Li+^He structure 
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Intrinsic structure of of '^^Be. 



Figure (a) shows a surface cut of the matter density 
p > 0.16 nucleons/fm'*. The density is calculated for the intrinsic AMD wave function before 
spin-parity projection. Density distributions for the single particle orbits of the valence neutrons 
are presented in the figures (b) and (c). 



in the states in tlie second K'^ = 3/22 band whicli starts from the J =3/22 state 
about 5 MeV higher than the ground state. The clusters of K'^ = 3/22 band are 
also exotic nuclei ^Li and ^He, and develop as remarkably as those in the K'^ = O2 
band of ^^Be. The mechanism of the ^He clustering development in ^"''B can be also 
described as the correlation among 4 neutrons in the sd-shell with 2 protons in the 
p-shell as similar way as -"^^Be. That is to say that these excited states with ^Li+^He 
clustering are based on the 4p — 2h state of the neutron p-shell. The spatial orbits 
of 4 valence neutrons are very similar as those of ^^Be presented in Figs 3(b) and 
3(c). The correlation among 4 neutrons with 2 protons is very unique even though 
the correlation among two neutrons and two protons has been sometimes discussed 
in relation with an a cluster. In this paper, we suggest one of the novel feature that 
4 neutrons in the sd-shell with 2 protons in the p-shell correlate to form a developed 
^He cluster in the neutron rich nuclei. It is an exciting problem whether or not the 
^He cluster formation may appear in other neutron-rich nuclei. 

Since the present results are obtained within the bound state approximation, 
particle decay width has not been discussed though they are important for the sta- 
bility of the excited states. Newly predicted states in K = O2 of ^^Be are expected 
to decay into the ^He-|-^He channel because of the developed clustering. We have 
calculated the partial decay width of ^^Be only for the *^He+®He channel with a 
method of reduced width amplitude, and found that the widths of 62 and 8^ states 
are less than 50 keV. For more detailed analysis of the total width including neutron 
decays and excited He decays, we require other microscopic frameworks such as a 
complex scaling method. 

§4. Shape coexistence of ^^Ar, ^''Ca 

In the previous section, structures of light neutron-rich have been discussed. 
Also in neutron-rich sd-shell nuclei we know many interesting phenomena such as 
vanishing of neutron magic number AT = 20 in ^^Mg. For study of the structures of 
neutron-rich nuclei, it is very useful to make systematic study associating unstable 
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nuclei with stable nuclei. For example, let us consider a problem of p-shell nuclei, the 
magic number N = 8 disappearance in ^^Be. The intruder ground state is considered 
to have 2 particles and 2 holes in neutron p-shell. We have already known the similar 
neutron configuration in the O2 state of ^^O which is described in terms of 4p — 4/i 
state with ^^C+a clustering. It is very interesting to imagine the structure change 
when the system varies from ^^0(02^) to ^■^Be with decrease of proton number of the 
system. In case of sd-shell nuclei, we think that it is helpful to compare the structures 
of stable nuclei near ^'^Ca with those of '^^Mg to solve the mechanism of neutron magic 
number vanishing in neutron-rich sd-shell nuclei. Fortunately many rotational bands 
of ^^Ar and '^'^Ca have been discovered with the recent 7-ray measurements. In this 
section, the study of the excited states of sd-shell nuclei '^^Ar and ^'^Ca is presented. 
We discuss the shape coexistence problem of these nuclei which is one of the recent 
hot subjects. 

For the structure study of sd-shell nuclei such as ^^Ar and ^'^Ca, we make vari- 
ation after parity projection but no spin projection (VBP) with AMD under a con- 
straint on the total oscillator quanta as mentioned in Sec.§. The interactions with a 
finite-range 3-body force are adopted because the ordinary effective forces such as the 
MVl force and the Gogny force are not appropriate to describe the binding energies 
and radii of nuclei covering wide mass number region from a to ^^Ca simultaneously. 
The central part of present interactions are as follows, 



Ventral = J:V^'^+ E ^^'^ (4) 

i<j i<j<k 



= (1 _ 




+ y,exp[-(^)2], (6) 

= VdeM-d{rl2 + 4s + 4,)% (7) 

Va = -198.34 MeV, Vb = 300.86 MeV, Vc = 22.5 MeV, (8) 

ra = 1.2 fm, n = 0.7 fm, = 0.9 fm, Va = 600 MeV, d = 0.8 fm'^ (9) 

m = 0.193,6 = -0.185, /i = 0.37, (10) 

where the ranges and strength parameters are chosen so as to reproduce reasonably 
the sizes of a and '^'^Ca, and the binding energies of a, ^^O and ^'^Ca. To choose 
interaction parameters the other important features like a+a phase shift and satu- 
ration property of the symmetric nuclear matter have been taken into consideration. 
In the total interaction, the spin orbit force of G3RS with the strength uis = 2500 
MeV and coulomb force are added to the central force. 

First the excited states of '^^ Ar are studied by VBP calculations with a constraint 
on the total oscillator quanta M. Here AM is defined by the deviation from the 
minimum oscillator quanta of the system; AM = M — Mmin, where Mmin is 52 in case 
of ^^Ar. After spin projection, we obtain an energy surface as a function of AM. It 
is found that there is a minimum point in the energies of 0"*" states at AM = 1, which 
corresponds to the ground state with a normal configuration. By diagonalyzing with 
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the ground state, we obtain the excited state with 4:hoj configurations at AJ\f = 4 
as a local minimum state. As a result, some rotational bands are constructed. In 
the ground K = band the intrinsic state has a normal oblate deformation. On 
the other hand the largely deformed state with 4fiuj configurations makes excited 
rotational bands K = 0^ and K = 2^ . In spite of the highly excited configurations 
the excited K = 0^ band starts from low energy region about 5 MeV and reaches high 
spin states J > 16. This excited K = band corresponds well to the experimental 
data observed recently with 7-ray transitions measurementsIlM'. The intrinsic state of 
the excited band deforms prolately as /3 = 0.3. It is interesting that another shape 
may coexist at AJ\f = 5. Namely, the triaxial shape with pentagon component 
appears to make another = 10'*' state at about 10 MeV in present calculations. 

In case of ^'^Ca, the calculations suggest that various kinds of shapes coexist 
in low energy region although '^''Ca is a doubly magic nucleus. First we find the 
ground state with Ohu configuration at AM = 2. It has spherical shape (a) as 
is expected. With the increase of oscillator quanta AM, we obtain energetically 
stable states with ifiuj configurations at AM = 5. One has an oblate shape (b) 
and the other has a prolate shape (c). Both intrinsic states construct rotational 
bands as shown in Fig. 4. In the figure only the excitation energies of the states 
with J < 16 are presented because of the limitation of computational calculations of 
spin projection. As the value AM increases, finally a super deformation (d) appears 
at AM = 10. By analyzing the single particle wave functions the super deformed 
state is found ot be dominated by 8huj configurations. The deformation parameter 
is ;9 ~ 0.4. It is surprising that the rotational band K = 0'^ given from the super 
deformation starts from the low energy region about 8 MeV excited in spite of the 
highly excited configurations of the doubly magic nucleus. As a result, many shapes 
coexist in ^'^Ca. It is suggested that the spherical ground state, the oblate excited 
state, the normal prolate excited state, and the prolate super deformed state appear 
in low energy region. Actually there exist many rotational bands in the experimental 
data measured recently The calculated super deformed band is considered to 
be the recently discovered rotational band from 5.273 MeV in the experimental 
data. Although the level spacing of the super deformed band is overestimated by 
present calculations, the results should be improved by extended calculations such 
as cranking methods and GCM calculations along the constraint AM. 

In the intrinsic states of the excited states, very strange shapes are found. Figure 
5 shows the density distributions of the intrinsic states of the excited bands. A 
hexagon like three leaves clover shape appears In the oblate state (b), while the 
super deformed state (d) has a parity asymmetric shape like a pear. One of the 
origins of such exotic shapes is clustering effects due to ^^C clusters. 

As mentioned above, many exotic shapes in the positive parity states are sug- 
gested in the present calculations. The appearance of a pear-like shape is an exciting 
prediction because the existence of parity asymmetric shapes in such heavy nuclei 
as ^"^Ca is an open problem. In order to confirm the parity asymmetric shapes it is 
very important to find parity doublet states in negative parity states. As is expected 
we find the parity doublet states of the super deformation with a pear shape. Also 
in the negative parity states, many shapes are found in the theoretical results. We 
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Fig. 4. Theoretical results of energy levels of positive parity 
states of ''"Ca. 



Fig. 5. Density distributions of the 
intrinsic states of the rotational 
bands. Figures (b), (c) and (d) 
correspond to the oblate shape, 
the prolate shape and the super 
deformation, respectively. 



consider the obtained K = 0~ at AM « 8 is the parity doublet of the super deformed 
band (d) because the intrinsic state is similar to each other. The excitation energy 
of the band head 1~ state is 12 MeV which is about 3 MeV higher than the energy 
of 0"*" state in the super deformation band. Present results indicate the possibility 
of the parity asymmetry shape due to clustering effects in heavier system. 

§5. Summary 

In summary, we have studied the structure of the excited states of ^■^Be, ^^Be 
and ^^B with variation after spin-parity projection based on the framework of AMD 
method. We have also reported our recent study of shape coexistence problems 
in ^^Ar and ^'^Ca with a constraint AMD. We have obtained a variety of exotic 
structures concerning with clustering. 

Various clustering structures in the excited states have been found in the theo- 
retical results in which we have succeeded to reproduce well the excitation energies 
of many levels of ^^Be. The breaking of the neuron magic number = 8 in ^^Be 
has been found. The recently measured excited states are described by the third 
K'^ = 0^ band which has a foreign clustering structure with developed ^He and ^He 
clusters. Exotic clustering structures such as ^He-|-®He and ^Li-|-^He have been pre- 
dicted in ^^Be and ^^B, respectively. The formation of the remarkable ^He clusters 
can be understood by an unique idea of the new type correlation among 4 neutrons 
in sd-shell and 2 protons in p-shell. This is the first full microscopic calculation 
which shows the exotic clustering in the excited states of ^^Be, ^"^Be and ^^B. The 
results suggest an important feature that the exotic clustering structures may exist 
very often in the excited states of neutron-rich nuclei. 

Structures of ^^Ar and ^'^Ca have been studied. In the results it is suggest that 
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many kinds of shapes coexist in these nuclei. In Ca, the spherical ground state, 
the normal prolate state, the oblate shape, and the prolately super deformation have 
been found. The prolate state and oblate state arc dominated by Ahuj configurations, 
while the super deformation originates from 8huj configurations. It is surprising 
that the rotational bands made of such highly excited configurations start at low 
excitation energies. The rotational band of the super deformation corresponds to 
the experimentally measured band which were observed recently in 7 ray transitions. 
In the theoretically obtained intrinsic states, very strange shapes are predicted in 
the excited states. It is very interesting that parity asymmetric shape like a pear 
is proposed in the super deformation. Although the negative parity bands have not 
been confirmed experimentally yet, present results predict that the parity doublet 
band with negative parity may exist at about 3 MeV higher than the positive parity 
band. 

If we vary the system from ^°Ca to ^^Mg by decreasing the proton number, how 
do the structures of the excited states of stable nuclei change ? Docs the excited 
state with 4p — Ah in the neutron shell appear in ^^Mg ? It is a very interesting 
subject to study structures of intruder states in neutron-rich nuclei associating with 
the excited of stable nuclei. 
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